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Abstract 

We present the results of a Bayesian analysis of a Regge model for K + A photoproduction. 
The model is based on the exchange of K + (494) and A"* + (892) trajectories in the t-channel. 
For different prior widths, we find decisive Bayesian evidence (A\n2 ~ 24) for a K + A photo- 
production Regge model with a positive vector coupling and a negative tensor coupling constant 
for the A"* + (892) trajectory, and a rotating phase factor for both trajectories. Using the \ 2 
minimization method, one could not draw this conclusion from the same dataset. 
Keywords: Regge phenomenology; Bayesian inference; Model selection 
PACS numbers: 11.55.Jy, 12.40.Nn, 13.60.Le 

1 Introduction 

The study of electromagnetic open strangeness or kaon-hyperon (KY) production is a key step 
towards understanding the structure of the nucleon[l, 2]. The focus of such studies has recently 
shifted from parameter estimation to model comparison, in part due to the importance of identifying 
the set of contributing resonances[3, 4, 5]. The statistical tools, however, have not been adapted 
to this new objective. The least-squares method in particular has often been stretched beyond its 
limits, being used not only as an optimization tool, but also as a model selection criterion. We 
advocate the Bayesian evidence as a more robust and well-founded tool for model comparison, and 
apply it to a Regge model for K + A photoproduction. 



2 Bayesian analysis 

The quantity we will use to quantify model fitness is the Bayesian evidence or marginal likelihood 
Z, which is proportional to the posterior probability P(M\ {dk}) of a model M given data {c4}[6]. 
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Z can be calculated by integrating over the model parameters olm 

Z = P({d k }\M) = f P({d k }\a M ,M) P{oc M \M) da M - (1) 

J V v ' V v ' 

C(a M ) 7r(a M ) 

We employ a uniform prior tt(cxm) and use the following limiting distribution to approximate the 
likelihood function [7] 

ri \ ~ 1 (x 2 ((*m) - k) 2 

c{aM) *^VTk exp it • (2) 

To compute the Bayesian evidence integrals, we employ the Nested Sampling (NS) method, de- 
veloped by Skilling[8]. The results can be interpreted qualitatively using Jeffreys' scale [9], i.e. 
AlnZ > 5 decisively favors the model with the highest evidence. 



3 Bayesian analysis of a Regge Model 

The Regge model under investigation is based on the exchange of the if + (494) and if*+(892) 
trajectories [10, 11]. The amplitude is derived from the i-channel Feynman amplitude by replacing 
the Feynman propagator by the respective Regge propagator, as described in Ref. [5]. The model's 
free parameters are the so-called sign factors (either 1, constant phase or e ~ l7ra (*- ) , rotating phase) 
and three continuous parameters, namely the strong coupling constant gK+A P 01 the K + trajectory 
and the tensor and vector couplings of the K* + trajectory, G^*» + = eg v ^, + Kp n K + K *+ /Air. 

Optimization of these parameters against the 72 available high-energy (E 7 > 5 GeV) data 
points[12, 13, 14] reveals that there are several model variants with comparable x 2 values [5]. The 
sign factors (e~ 47rQ (*) or 1) as well as the signs of G V K , + and G l K *+ cannot be established conclusively 
using the x 2 " mc thod. These sign and phase ambiguities may not seem important for the Regge 
model itself. However, for the Regge- plus- resonance (RPR) model, in which the Regge background 
is complemented with s-channel nucleon (N*) resonances, an exact determination of the background 
parameters is of major importance as it affects the extraction of the resonance information. 



4 Results 

We employ a uniform prior for all continuous parameters and allow for a 20% deviation from the 
SU(3) prediction for the coupling constant gx+A P when determining its prior interval[15]. For 
G^, + , we compare results using prior widths of 100, 1000 and 10000 to ensure that the results 
do not depend on the choice of the prior. Indeed, multiplying the prior width with a factor of 
10 results in a difference of less than 5% in the computed values of InZ. Thereby, one should 
bear in mind that a reduced sampling efficiency leads to considerably larger errors. Repeating the 
calculations assuming 40% SU(3) symmetry breaking does not significantly influence the results. 
More importantly, the ranking of the models is not affected by the above prior modifications [16]. 

We find that the value of InZ for the best model variant for K + A photoproduction, with a 
rotating phase for both K + and K* + trajectories, positive vector and negative tensor coupling, 
is 24.30 ± 0.75 above the second-best variant. This result decisively resolves the sign and phase 
ambiguity for K + A photoproduction, which could not be achieved in a previous fit of the model to 
the high-energy data using the x 2 method [5]. 



5 Conclusions and outlook 

Bayesian inference provides us with an excellent tool for model comparison. We have demonstrated 
this by using the Nested Sampling algorithm to compute the Bayesian evidence for different model 
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variants of a Regge model for K + A photoproduction. The Nested Sampling method has many ap- 
plications, both for the RPR model and for other research. One of these applications is the accurate 
estimation of model parameters as well as the elimination of nuisance parameters. More impor- 
tantly, however, this method may provide us with a means to address the missing-resonance problem 
by calculating the probability of individual resonance contributions in a Bayesian framework. This 
is an approach we intend to explore in the near future. 
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